Fractal Conductance Fluctuations of Classical Origin 
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In mesoscopic systems conductance fluctuations are a sensitive probe of electron dynamics and 
chaotic phenomena. We show that the conductance of a purely classical chaotic system with either 
fully chaotic or mixed phase space generically exhibits fractal conductance fluctuations unrelated 
to quantum interference. This might explain the unexpected dependence of the fractal dimension 
of the conductance curves on the (quantum) phase breaking length observed in experiments on 
semiconductor quantum dots. 
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A prominent feature of transport in mesoscopic sys- 
tems is that the conductance as a function of an external 
parameter (e.g. a gate voltage or a magnetic field) shows 
reproducible fluctuations caused by quantum interference 
PJ . A prediction from semiclassical theory that inspired a 
number of both theoretical and experimental works in the 
fields of mesoscopic systems and quantum chaos was that 
in chaotic systems with a mixed phase space these fluc- 
tuations would result in fractal conductance curves 0,0]. 
Such fractal conductance fluctuations (FCFs) have since 
been confirmed in gold nanowires and in mesoscopic semi- 
conductor billiards in various experiments 0, H, 0, 0] ■ 
In addition FCFs have more recently been predicted to 
occur in strongly dynamically localized [9( and in diffu- 
sive systems ■ Due to the quantum nature of the FCFs 
it came as a surprise when recent experiments indicated 
that decoherence does not destroy the fractal nature of 
the conductance curve but only changes its fractal di- 
mension [ill, [12] ■ In the present letter we show, that 
the conductance of purely classical (i.e. incoherent) low- 
dimensional Hamiltonian systems very fundamentally ex- 
hibits fractal fluctuations, as long as transport is at least 
partially conducted by chaotic dynamics. Thus mixed 
phase space systems and fully chaotic systems alike gen- 
erally show FCFs with a fractal dimension that is de- 
termined analytically. We show that it is governed by 
fundamental properties of chaotic dynamics. 

In a disordered mesoscopic conductor - which is 
smaller than the phase coherence length of the charge 
carriers but large compared to the average impurity spac- 
ing - the transmission is the result of the interference of 
many different, multiply-scattered and complicated paths 
through the system. As these paths are typically very 
long compared to the wave length of the charge carriers, 
the accumulated phase along a path changes basically 
randomly when an external parameter such as the en- 
ergy or the magnetic field is varied. This results in a ran- 
dom interference pattern, i.e. reproducible fluctuations 
in the conductance of a universal magnitude on the order 
of 2e 2 /h, the so called universal conductance fluctuations 
(UCFs). For a review see (T^l or The role of disorder 
in providing a distribution of random phases can as well 



be taken by chaos. Thus ballistic mesoscopic cavities like 
quantum dots in high mobility two-dimensional electron 
gases that form chaotic billiards show the same universal 
fluctuations [3, El E3]- If the average of the phase gain 
accumulated on the different paths traversing the system 
exists, the conductance curves are smooth on parameter 
scales that correspond to a change of the average phase 
gain on the order of and smaller than the wave length of 
the carriers. In systems with mixed phase space, where 
chaotic and regular motion coexist, this phase gain, how- 
ever, is typically algebraically distributed and its aver- 
age phase gain does not exist (neglecting the finiteness 
of the coherence length and assuming the semiclassical 
limit H e ff — > ; the role of the finite H e ff is discussed in 
H), Therefore, as shown in [§], the conductance curve 
of such a system fluctuates on all parameter scales and 
forms a fractal. The fractal dimension D is connected 
to the exponent 7 of the algebraic distribution of phase 
gains by D = 2 — ^. 

Experiments on quantum dots that study the depen- 
dence of the conductance fluctuations on several system 
parameters like size and temperature seem to partly 
tradict the semiclassical theory of fractal scaling 
Namely it was found that with decreasing coherence 
length the scaling region over which the fractal was ob- 
served did not shrink - as would be expected from the 
semiclassical arguments -, but that the fractal dimen- 
sion changed. An implicit assumption of the semiclas- 
sical theory is that the classical dynamics remains un- 
changed as the external parameter is varied and thus 
only phase changes are relevant. In most experimental 
setups, however, the classical phase space changes with 
variation of the control parameter. In this article we show 
that the classical chaotic dynamics itself already leads to 
fractal conductance curves! Moreover, from this follows 
that even on very small parameter scales the fluctuations 
due to changes in the classical dynamics are important. 
In general the conductance curve is a superposition of 
two fractals: one originating in interference which is sup- 
pressed by decoherence to reveal the fractal fluctuations 
reflecting the changes in the classical phase space struc- 
ture. In addition, we predict that FCFs are not only 
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observable in systems with a mixed phase space but in 
purely chaotic systems. 
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Figure 1: Classical Conductance g{B) through a stadium 
(left, geometry as in ref. 0]) and square billiard (right, ge- 
ometry as in ref. [HI]) versus magnetic field B. Both fluc- 
tuating conductance curves are fractals (see insets and text). 
Their respective dimensions are D « 1.28 for the stadium and 
D ~ 1.25 for the square billiard. The fractal dimensions are 
in good agreement with experimental measurements 0, EH- 

As a starting point of our investigations and to connect 
it to the experiments we numerically study the classical 
conductance through a rectangle (hard-wall) and a sta- 
dium billiard (soft-wall) as a function of a magnetic field 
as shown in Fig. [TJ (Throughout this article, we will 
study the transmission, which, in accordance with the 
Landauer theory of conductance, is proportional to the 
conductance, see e.g. [13].) Note that not only the phase 
space of the stadium but also of the rectangle billiard is 
mixed in the presence of a perpendicular magnetic field. 
In both modified version (Til . of the box- 

counting analysis clearly reveals the fractal nature of the 
conductance curves. As the simulation is purely classi- 
cal, the fractal scaling cannot be caused by interference 
effects. So what is the underlying mechanism for the 
fractality of the conductance curve and how can we un- 
derstand its dimension? 

To study this mechanism in detail we will, because 
of its numerical and conceptual advantages, analyze the 
transport in Chirikov's standard map 

HElIIEi. This 

paradigmatic system shows all the richness of Hamilto- 
nian chaos. And since - as will become apparent below - 
our theory relies only on very fundamental properties of 
chaotic systems, it is a natural choice as a model system. 
The standard map is defined by 

9' = p + d 

p' = p + K sin 9' 

with momentum p , angle 9 and the 'nonlinearity pa- 
rameter' K, which drives the dynamics from fully inte- 
grable (K = 0) to fully chaotic (K > 7), In between 
the phase space is mixed. The standard map can be seen 
as the Poincare surface of a conservative system of two 
degrees of freedom. As such the map can by viewed to di- 
rectly correspond to the Poincare map at the boundary 
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Figure 2: How lobes translate into fluctuations: In the lower 
row the entryset of the standard map with absorbing bound- 
ary conditions at ±37r for K = 7.5 and K = 7.6 resp. can be 
seen. The three pictures in the center row show the magnifi- 
cation of the central sections of the entryset for three different 
values of K = 7J>, 7.55 and 7.6. The transmission T(K) for 
K = 7.5 . . . 7.6 [20| | is shown in the top left picture. Note that 
a small change in K shifts the lobes vertically, but conserves 
the overall phase space structure, and that the largest fluc- 
tuations are caused by intersection with the apex of lobes. 
Starting from K — 7.5, a large transmission lobe is cut by 
the horizontal line (see text), i.e. the transmission increases 
with K. In the same way, e.g. the fluctuations of T(K) near 
K — 7.55 can be understood. The box-counting analysis re- 
veals a fractal structure (top right). 



of a chaotic ballistic cavity, connecting it conceptually 
with the experimental system. We introduce absorbing 
boundary conditions (see e.g. ref. [3]), be. when p ex- 
ceeds (drops below) a maximum (minimum) threshold 
value, the particle is transmitted (reflected) and leaves 
the cavity. As can be seen right from the definition of 
the standard map, the envelope of the entryset (which is, 
the phase space projection of the injection lead) is simply 
half a period of a sine function times K. 

A trajectory entering the system eventually con- 
tributes either to the total transmission or reflection, and 
we mark the corresponding point in the entryset by a 
color code (transmission: red, reflection: blue). Chaotic 
dynamics, through its fundamental property of stretch- 
ing and folding in phase space, leads to a lobe structure 
(see Fig. [2] (bottom)), which is typical for chaotic systems 
and not special to the standard map. The distribution 
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of widths w of lobes exhibits a power law 

n(w) oc w~ a . 



(1) 



The lobe structure is translated into transmission by 
summing up the intersections of the transmission lobes 
along a horizontal line, see Fig. [2 A lobe of thickness w 
gives rise to a maximum contribution AT oc w@ . Varia- 
tion of the external parameter K leads to a fractal trans- 
mission curve T(K) with D w 1.25. 

How does the fractal dimension depend on the power 
law distribution of lobe-widths and the curvature of the 
lobes? To this aim, we study a random sequence of curve 
segments mimicking the intersection of consecutive lobes 
of widths w, distributed algebraically with exponent a 
and curved like w 13 . We define Xi := Y?j=i w j an d 



T{X) = {-l)\X-X l ) 



Xi < X < X, 



An example of this curve of "random lobes" with a = 1.9 
and /? = g is shown in Fig. [3] (top) . The box-counting 
analysis clearly reveals a fractal structure. 
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Figure 3: Transmission T(X) for lobes (red upper curve, 
shifted along the y-axis for clarity) and stripes (black lower 
curve) for one and the same random distribution with a = 1.9, 
ft = 0.5. The inset shows the box-counting analysis for the 
upper (red triangles) and lower transmission curve (black 
squares). The regression line is drawn for the upper curve, 
whose fractal dimension is 1.41. 

We further simplify the problem by replacing the lobes 
by a sequence of stripes of widths x with power law dis- 
tribution n(x) oc x a . Dispensing with the sign of the 
fluctuation, the transmission reads 

T{X) = {X l+1 -Xif. 

This yields histogrammatic transmission curves T(X) 
like the bottom curve of Fig. 03 As shown in the in- 
set, the fractal dimension of the resulting transmission 
curve remains unchanged compared to the correspond- 
ing calculation with random lobes within the precision of 
the box-counting analysis. Thus, the fractal dimension 



of the curve does not change noticeably when considering 
stripes instead of lobes and also when neglecting the sign 
of each contribution, confirming the intuition, that the 
fractal dimension depends only on the relative scaling, 
i.e. a and (3, but not on the detailed form of the curve 
sections. 

For these curves like the bottom one of Fig. [3] with 
a — /3 > 1, we can give an analytical expression for the 
fractal dimension and then estimate the fractal dimension 
of the transmission curve in the standard map. We ap- 
ply the box-counting method, which we therefore review 
shortly (see e.g. (l9l | for a more detailed introduction). In 
this approach the fractal curve lying in a n— dimensional 
space is covered by a n— dimensional grid. Let the grid 
consist of boxes of length scale s. The box-counting di- 
mension is then given by 
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log N(s) 
s^o log(s) 



lim 



(2) 



where N(s) is the number of non-empty boxes. For 
our problem, we divide N(s) into three contributions 
N(s) = n a + rib + n c , as schematically drawn in Fig.[D[A). 
The number n a of vertically placed boxes (see mark (a)) 
covering contributions from stripes of widths x > s, reads 



1 f°° 

n a (s) oc - / p(x)x l3 dx oc s~ (Q-/3 ). 



(3) 



Secondly, the number rib of horizontally placed boxes 
covering horizontal contributions of stripes of widths 
larger than s, see Fig. H]A.(b), is given by 



n b (s) = - 



1 f°° 1 f°° 

- I p(x)xdx < — p(x)xdx. 



(4) 



Hence rib scales like s _1 and can be neglected in com- 
parison to n a because of a — (3 > 1. Finally, we de- 
termine an upper estimate for the number n c of verti- 
cally placed boxes covering the contribution from stripes 
of widths x < s. The total length of these widths is 
L(s) — J Q p{x)xdx, therefor L(s)/s boxes are needed to 
cover the length. Inflating all heights of the stripes x < s 
to the maximum possible size s 13 , see Fig. H]A.(c), we find 



. . L(S) _ a + g 

n c (s) < — — — oc s a+p . 

s s 



(5) 



For s<l thus the dominant terms is n a (s) . With Eq. ^ , 
N(s) gives rise to the box-counting dimension [25J 
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lo gs -°+l 
lim — = a — p. 

s^O log S 



(6) 



To connect the analytical result with the calculations 
of the transmission of the open standard map, we nu- 
merically estimate the distribution of lobe-widths in the 
entryset for K — 8, finding a w 1.9, as shown in Fig.[4j3. 
Together with (3 = | , corresponding to first order Taylor 
expansion of the cosine function, Eq. ((6]) predicts a fractal 
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Figure 4: A. Schematic transmission according to Fig. [3] (bot- 
tom), covered with boxes of size a. There are three contribu- 
tions marked (a-c). B. Total number Ni n t(w) = n(w')dw' 
of lobes (for the open standard map with \p\ < 4tt) of 
width larger than w on a double logarithmic scale. The 
four curves show estimates for increasing resolution 
10 -5 . . . 10 -8 . The curves clearly approach a power law corre- 
sponding to n(w) oc w~ 19 . The insets show the transmission 
curve T(K) for values K = 8.0 . . . 8.1 calculated from 2 x 10 13 
trajectories and its box-counting dimension. 

dimension D w 1.4. Direct analysis of the transmission 
curve (see insets of Fig. [4)3) yields a fractal dimension 



D sa 1.39, in good agreement with the expected value. 

How can a power law distribution of lobe widths 
emerge in a fully chaotic open system? One might rather 
expect to find an exponential distribution of lobes in a 
fully chaotic system. To see why the distribution is al- 
gebraic, however, let us examine the simplest case of an 
open chaotic area preserving map the dynamics of which 
is governed by a single, positive homogeneous Lyapunov 
exponent A. In each iteration phase space structures are 
stretched in one direction by exp(A), shrunk by exp(— A) 
in the other and then folded back. The entryset of the 
open system is thus stretched into lobes of decaying width 
w(U) oc exp(— XU). The phase space volume flux out 
of the system decays exponentially as it is typical for a 
fully chaotic phase space, i.e. T(ti) oc exp(— U/t), with 
(mean) dwelltime r. The area T(ti)At is the fraction 
of the exitset that leaves the system at time With 
ti(w) oc — ln( w)/A the number of lobes of width w in the 
exitset is [26| 

T(t l (w))At 1 M(w). i , 
#w = K y " — oc - exp(— = vj — - 1 . 
w w At 

This suggests that the power law distribution of lobe 
widths is a generic property even for fully chaotic sys- 
tems. A quantitative expression for the exponent, how- 
ever, is not as easy to derive, as e.g. the Lyapunov expo- 
nent for the standard map is not homogeneous. 

In conclusion, we have shown that transport through 
chaotic systems due to the typical lobe structure of 
the phase space in general produces fractal conductance 
curves, where the fractal dimension reflects the distribu- 
tion of lobes in the exit- /entryset. In contrast to the 
semiclassical effect the size of the fluctuations is not uni- 
versal but depends on specific system parameters. Due 
to the fractal nature of the classical conductance, how- 
ever, there is no parameter scale that separates coherent 
and incoherent fluctuations. 
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